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Abstract 

We study the observational constraints on the exponential gravity model of f(R) = —/3R S (1 — 
e~ R l Rs ). We use the latest observational data including Supernova Cosmology Project (SCP) 
Union2 compilation, Two-Degree Field Galaxy Redshift Survey (2dFGRS), Sloan Digital Sky 
Survey Data Release 7 (SDSS DR7) and Seven- Year Wilkinson Microwave Anisotropy Probe 
(WMAP7) in our analysis. From these observations, we obtain a lower bound on the model 
parameter f3 at 1.27 (95% CL) but no appreciable upper bound. The constraint on the present 
matter density parameter is 0.245 < ft^ < 0.311 (95% CL). We also find out the best-fit value of 
model parameters on several cases. 
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I. INTRODUCTION 



Cosmic observations from type la supernovae (SNe la) l|, |2j, large scale structure (LSS) 
3, 4], baryon acoustic oscillations (BAO) {|3| and cosmic microwave background (CMB) 
6|, LZj indicate that our universe is undergoing an accelerating expansion. The reason for this 
acceleration, the so-called dark energy problem, remains a fascinating question today. The 
simplest model to explain this problem is the ACDM model, in which a time independent 
energy density is added to the universe. However, the ACDM model suffers from both 



fine-tuning and coincidence problems |8Hl3|. In general, the ways to understand the cosmic 
acceleration can be separated into two branches. One is to modify the matter by introducing 
some kind of "dark energy" . The other one is to modify Einstein's general relativity - the 
modification of gravity. 

In modified gravity, one of the popular approaches is to promote the Ricci scalar R in 
the Einstein-Hibert action to a function, f(R). Although there are several viable f(R) 
models, many of them are restricted to the regimes to be effectively identical to the ACDM 
by the observational constraints. Recently, Linder [l^] has explored an f(R) theory named 



"exponential gravity", which has also been discussed in Refs. |15l-ll7l|. The exponential 
gravity has the feature that it allows the relaxation of fine-tuning and it has only one 
more parameter than the ACDM model. In addition, the exponential gravity satisfies all 



conditions for the viability 18] such as the local gravity constraint, stability of the late- 
time de Sitter point, constraints from the violation of the equivalence principle, stability of 
cosmological perturbations, positivity of the effective gravitational coupling, and asymptotic 
behavior to the ACDM model in the high curvature regime. In this paper, we will study 
the constraints given by latest observational data, reexamine the alleviation of the fine- 
tuning problem, and find the possibility of the derivation from ACDM. We use units of 
&b = c = h — 1 and the gravitational constant is given by G = M^ 2 with the Planck mass 
of M P1 = 1.2 x 10 19 GeV. 

The paper is organized as follows. In Sec. HH we review equations of motion and the 
asymptotic behavior at the high redshift regime in the exponential gravity model. In Sec. 
IHIt we discuss the observations and methods. We show our results in Sec. [TV] Finally, 
conclusions are given in Sec. [V] 



II. EXPONENTIAL GRAVITY 



The action of f(R) gravity with matter is given by 



= i / d^V 2 



2k 2 

where k 2 = 87rG and f(R) is a function of the Ricci scalar curvature R. In this paper, we 
focus on the exponential gravity model [ljj], given by 

f(R) = -PR s (l-e- R ' R °), (2.2) 

where R s is related to the characteristic curvature modification scale. Since the product of 
(3 and R s can be determined by the present matter density [ijj], we can choose (3 and 
Q^n as the free parameters in the model. 

We use the standard metric formalism. From the action (12.11) . the modified Friedmann 
equation of motion becomes [ljj] 

H 2 = ^ + UfRR -f)~H 2 (f R + f RR R>) , (2.3) 
3 o 

where H = a/a is the Hubble parameter, a subscript R denotes the derivative with respect 
to R, a prime represents d/dlna, and pu = pm + Pr is the energy density of all perfect 
fluids of generic matter including (non-relativistic) matter, denoted by m, and relativistic 
particles, denoted by r. Here, we only consider the matter density. Since the modification by 
the exponential gravity only happens at the low redshift, the contributions from relativistic 
particles are negligible. In a flat spacetime, the Ricci scalar is given by 

R = 12H 2 + 6HH'. 

Following Hu and Sawicki's parameterization [2(3], we define 

_ p DE H 2 _ 3 R _ 3 

Vh = 7. = — ? - a , y R = — - 3a , (2.4) 
p^ m A m A 

where m 2 = K 2 p° m /3, pde is the effective dark energy density, and p Q m is the present matter 
density. Then, Eqs. (12. 3 p and (12 ,4p can be rewritten as two coupled differential equations, 

Vh = f - 4y fl (2.5) 

and 
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H 2 



R 



+ 



/ 



6m 2 7 6m 2 



H 2 f RR 

where R and H 2 can be further replaced by y R and x/h from equations in (j22 
Eqs. (12. 5 p and ( 12. 6p . we obtain a second order differential equation of yn, 



Vh + JiVh + + J 3 = 0, 



(2.6) 
Combining 

(2.7) 



where 



Ji = 

Vh -f 
J 3 = -3a" 3 



1 



R 



y H + a 3 6m 2 f RR ' 
, /*a- 3 + //3m 2 



2/// + a' 



6m 2 / 



with 



i? = m 2 [3 (y' H + 4y H ) + 3a" 3 ] . 



(2- 



(2.9) 



Solving Eq. (12.71) numerically, we can get the evolution of the Hubble parameter in the low 
redshift regime (z = ~ 4). The effective dark energy equation of state wde is given by 



w DE = -1 - 



(2.10) 



In the high redshift regime (z > 4), the exponential factor e~ R / Rs of /(i?) in Eq. (12. 2p 
becomes negligible (e _iJ//ii:s < 10 -5 ). The exponential gravity model behaves essentially like 
a cosmological constant model with the dark energy density parameter Q\ = @Rs/6Hq = 
^mUHizhigh)- Thus, the Hubble parameter as a function of z in this regime can be expressed 
as 



H(z) = H Q JW (1 + zY + W (1 + zf + 



0Rs 
6tf 2 ' 



(2.11) 



where is the density parameter of relativistic particles including photons and neutrinos^. 
The equation (I2.1ip will be used in the data fitting of CMB and the high redshift part of 
BAO in section HTT1 



= 0° (1 + 0.2271iV e //), where fl° is the 



the effective number of neutrino species 



2l|. 



present fractional photon energy density and N e ff = 3.04 is 
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III. OBSERVATIONAL CONSTRAINTS 



To constrain the free parameters of f3 and in the exponential gravity model, we use 
three kinds of the observational data including SNe la, BAO and CMB. The SNe la and 
CMB data lead to constraints at the low and high redshift regimes, respectively, while the 
BAO data provide constraints at the both regimes. 



A. Type la Supernovae (SNe la) 

The observations of SNe la, known as "standard candles", give us the information about 
the luminosity distance asa function of the redshift z. The distance modulus /i is defined 

as 

fJ>th{zi) = 5)og 10 D L (zi) + /i , (3.1) 

where /io = 42.38 — 51og 10 h with Hq = h ■ 100km/ s/Mpc is the present value of the Hubble 
parameter. The Hubble-free luminosity distance for the flat universe is 

d -m = <- i+ (3 ' 2) 

where E(z) = H(z)/Hq. The \ 2 °f the SNe la data is 

2 _ \- [Vobs(Zj) - Hth(zj)] 2 
X-SN — / . ~2 ' \ 6 - 6 ) 

*~r erf 

where /i ft s is the observed value of the distance modulus. Since the absolute magnitude 
of SNe la is unknown, we should minimize y| y with respect to /io, which relates to the 



absolute magnitude, and expand it to be 22|, |23| 



Xj SN = A-2(M>B + ii*C t (3.4) 

where 

[Hobs(zi) - fJ*h(zi; /io = 0)] 2 



The minimum of Xsn with respect to /to is 

B 2 



X 2 sn = A- (3.6) 
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We adopt this Xsn f° r our later \ 2 minimization. We will use the data from the Supernova 
Cosmology Project (SCP) Union2 compilation, which contains 557 supernovae 24J, ranging 
from z = 0.015 to z = 1.4. 



B. Baryon Acoustic Oscillations (BAO) 



The observation of BAO measures the distance ratios of d z = r a (zn)/ 'Dy(z), where Dy 
is the volume-averaged distance, r s is the comoving sound horizon and Zd is the redshift at 
the drag epoch [25J. The volume- averaged distance Dy(z) is defined as [5] 

1/3 



Dv(z) 



[l + z) 2 Dl(z) 



Hz) 



(3.7) 



where Da(z) is the proper angular diameter distance: 

1 f z dz' 



D A z) 



l + zj H(z 



— , (for flat universe). 



(3.{ 



The comoving sound horizon r s (z) is given by 

I ri/(i+«) 
r s (z) = —= 



da 



V^o a?H{~' 



(3.9) 



where and Q° are the present values of baryon and photon density parameters, respec- 
tively. We use n° b = 0.022765/i- 2 and = 2.469 x lO" 5 /*" 2 (21]. The fitting formula for z d 
is given by 26] 

1291(^/i 2 ) - 251 



Zd 



1 + 0.659(^/i 2 ) - 828 



[l + &i(^ 2 ) 62 ] 



(3.10) 



where 



61 = 0.313(^/i 



z,2\-0.419 



[1 + 0.607(^/ i 2 ) - 674 ] , 



0.238(fiyi 



l2\0.223 



(3.11) 



The typical value of z d is about 1021 with = 0.276 and h = 0.705. Since z d is in the 
high redshift regime, we use Eq. (12. lip to calculate r s {zd). On the other hand, Dy(z) is 
evaluated by the numerical result of Eq. ( 12. 7p as it is in the low redshift regime. 

The BAO data from the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and the 
Sloan Digital Sky Survey Data Release 7 (SDSS DR7) [25] measured the distance ratio d z at 
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two redshifts z = 0.2 and z = 0.35 to be d°% 2 = 0.1905±0.0061 and d° z % 35 = 0.1097±0.0036 
with the inverse covariance matrix: 

/ 30124 -17227 \ 
C B \o = • (3-12) 

\ -17227 86977 J 

The x 2 f° r the BAO data is 

Xbao = ( x i^BAo ~ x i b BAo)(C BAO )ij(x t J ^ BAO — x° b BAO ), (3.13) 
where x i)B Ao = (4.2, ^0.35)- 



C. Cosmic Microwave Background (CMB) 



The CMB is sensitive to the distance to the decoupling epoch 27J. It can give con- 
straints on the model in the high redshift regime (z ~ 1000). The CMB data are taken from 
Wilkinson Microwave Anisotropy Probe (WMAP) observations 2l|. To use the WMAP 

data, we compare three quantities: (i) the acoustic scale I a, 

7tD a (z*) 



Ia(z*) = (1 + 2*)- 



(ii) the shift parameter R 



^H (l + z*)D A (z*), 



(3.14) 



(3.15) 



and (iii) the redshift of the decoupling epoch z*. The fitting function of z* is given by 29] 

(3.16) 

where 



fji 



1048 [1 + 0.00124(n°fr 2 )- a738 ] [1 + gi {Q° m h 2 ) 92 ] 

0.0783(iy)-°- 238 0.560 
= l + 39.5(fi°/i 2 ) - 763 ' 



92 



l + 21.1(fig/i 2 ) 1 - 81 ' 



The x 2 of the CMB data is 



XCMB 



X 



til 



„obs 



CMB X j,CMB)i 



(3.17) 



(3.18) 



— 'r obs Vr -1 (rr th 
,CMB • L i,CMB)V^CMB)'i-3\- L 3.(. 

where x^cmb = (Ia(z*), R(z*), z *) and C^ l MB is the inverse covariance matrix. The data 
from Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP 7) observations 2lJ lead 
to Ia(z*) = 302.09, R(z*) = 1.725 and z* = 1091.3 with the inverse covariance matrix: 

I 2.305 29.698 -1.333 ^ 

29.698 6825.27 -113.180 . (3.19) 
-1.333 -113.180 3.414 



^ CMB 



\ 



) 



Finally, the \ 2 °f an t ne observational data is 

X 2 = X 2 sn + xIao + Xomb- (3- 20 ) 

In our fitting process, we did not use the Markov chain Monte Carlo (MCMC) approach 
because the numerical calculation for each solution of f (R) theory is very time-consuming, 



and the necessary change to the code like CosmoMC |30j is very extensive with no obvious 
benefit in our study of the exponential gravity. Therefore, we take the simple \ 2 method 
as our main fitting procedure. The ACDM result obtained from SNe la, BAO and CMB 



constraints with this \ 2 method is 0,^ = 0.276^^013) whi 



3l|. We note that the fitting in Ref. 



e that with the MCMC method is 



3 II ] has also included the observa- 



OO _ n O79+0.013 
"m — u - z ' z -0.011 

tional constraints from the radial BAO and Hubble parameter H(z). 

TABLE I. The best-fit values of the matter density parameter (68% CL) and x 2 f° r the 
exponential gravity model with j3 = 2,3,4 and the ACDM model. Note that the error for fJ^ is 
obtained when f3 is fixed. 



Model 






0° 


x 2 




P = 


= 2 


2?4 +0.014 
u - z '^0.013 


546.7136 


Exponential Gravity 


P = 


= 3 


276+ ' 014 


545.3836 




P = 


= 4 


276+ - 014 
u.z(u_ 013 


545.1721 


ACDM 








545.1522 



IV. RESULTS 

Based on the methods described in Sec. IHH we now examine the parameter space of 
the exponential gravity model. In Fig. [TJ we present likelihood contour plots at 68.3, 95.4 
and 99.7% confidence levels obtained from the SNe la, BAO and CMB constraints. The 
results show that the observational data give no upper bound on the model parameter (3, 
making it a free parameter. Hence, there is no fine-tuning problem. However, a larger value 
of /3, which is closer to the ACDM model, is slightly preferred by the observational data 
as expected. The lower bound on f3 is /3 > 1.27 (95% CL). The present matter density 
parameter is constrained to 0.245 < < 0.311 (95% CL), which agrees with the 
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FIG. 1. The 68.3%, 95.4% and 99.7% confidence intervals for the exponential gravity model, 
constrained by the SNe la, BAO, and CMB data. The best-fit point in this parameter region is 
marked with a plus sign. 

current observations. The best-fit value (smallest x 2 ) m the parameter space between /3 = 1 
and cl is x 2 = 545.1721 with (3 = 4 and fi^ = 0.276. The comparison of the best-fit tt° m 
and x 2 f° r the model with — 2,3,4 and ACDM is shown in Table [H 



2 We only concentrate on the region of 1 < /3 < 4. For j3 > 4, it is almost the ACDM model. For /? < 1, it 
is ruled out by the local gravity constraints and the stability of the de-Sitter phase. 
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In Fig. |2l we illustrate the evolution of the effective dark energy equation of state wde for 
(3 = 2,3,4 with their best-fit which is given in Table [H We can see that, for every value 
of /3, the effective dark energy equation of state woe starts at the phase of a cosmo logical 
constant Wde — — 1 and evolves from the phantom phase (w^e < — 1) to the non-phantom 
phase (woe > — 1)- And, for larger value of /3, the deviation from cosmological constant 
phase (wde = —1) become smaller. For /3 — 2, there is still another small oscillation after 
the main phantom phase crossing. Negative z means the future evolution. It is clear that 
the exponential gravity model has the feature of crossing the phantom phase in the past as 



well as the future 



32| 



In Fig. |3l we depict the effective dark energy density Qde and non-relativistic matter 
density Q m vs. the redshift z. 
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FIG. 2. Evolution of the effective dark energy equation of state wde corresponding to (3 = 2,3,4 
with their best-fit $7^ given in Table HI 
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z 



FIG. 3. The evolutions of the effective dark energy density parameter VLde and non-relativistic 
matter density parameter Vt m as functions of z, where the solid lines indicate the exponential 
gravity model with f3 = 1.27 and the best-fit = 0.270 and the dashed lines represent the 
ACDM model with Qj^ = 0.276. For a higher value of (3, the evolution becomes closer to that in 
ACDM. 

V. CONCLUSION 

We have studied the exponential gravity model. In the low redshift regime, we follow Hu 
and Sawicki's parameterization to form the differential equation for the exponential gravity 
and solve it numerically. In the high redshift regime, we take advantage of the asymptotic 
behavior of the exponential gravity toward an effective cosmo logical constant. The analytical 
form of the Hubble parameter if as a function of the redshift z can be expressed in the high 
redshift limit. We have constrained the parameter space of the model by the SNe la, BAO 
and CMB data. We have found that there is a lower bound on the model parameter f3 
at 1.27 but no upper limit, and fl^ is constrained to the concordance value. This means 
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that the exponential gravity model shows no need of fine-tuning. Nevertheless, the ACDM 
model is still included by the observational constraints since /3 — > oo corresponds to the 
model. Current observational data still lack the ability to distinguish between the ACDM 
and exponential gravity models. 

Finally, we remark that as seen from Fig. El the noticeable difference between the ex- 
ponential gravity and ACDM models lies in the regime 0.2 < z < 1, and is maximized 
at z — 0.5 if we compare their expected distance modulus. An improvement on the BAO 
observation may give a stronger constraint on this redshift regime or higher. The ongoing: 
and future dark energy survey projects which will observe BAO include WiggleZ 



33[, BOSS 



(Baryon Osci 



lation Spectroscopic Survay) [34J, HETDEX (Hobby- Eberly Dark Energy Ex- 



periment) 35|, EUCLID 36_|, JDEM (Joint Dark Energy Mission) /Omega with Wide Field 



Infrared Survey Telescope (WFIRST) js^j, BigBOSS (Big Baryon Oscillation Spec-troscopic 



Survay) |38j, SKA (Square Kilometer Array) [39|, LSST (Large Synoptic Survey Telescope) 



40| and DES (Dark Energy Survey) j4l|. In addition, it is known that the measurement on 



the growth rate of f g (z) = dhx8 m /dlna has the potential to distinguish the models with the 
same expansion history but different physics. In the exponential gravity case, the growth 
index is 7 = 0.540 for /3 = 2. It is clear that if those surveys such as WiggleZ, EUCLID, 
BigBOSS and JDEM/Omega can measure the growth rate with a high accuracy, they will 
be able to discriminate the exponential gravity from the ACDM model. 

ACKNOWLEDGMENTS 



We thank Dr. K. Bamba for many helpful discussions and suggestions. The work is sup- 
ported in part by the National Science Council of R.O.C. under: Grant NSC-98-2112-M- 
007-008-MY3 and National Tsing Hua University under the Boost Program #: 97N2309F1. 



[1] A. G. Riess et al. [SNST Collaboration], Astron. J. 116, 1009 (1998). 

[2] S. Perlmutter et al. [SNCP Collaboration], Astrophys. J. 517, 565 (1999). 

[3] M. Tegmark et al, [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004). 

[4] U. Seljak et al. [SDSS Collaboration], Phys. Rev. D 71, 103515 (2005). 

[5] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005). 



12 



D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003). 

D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007). 
S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 

V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000). 
S. M. Carroll, Living Reviews in Relativity 4 (2001). 
P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003). 
T. Padmanabhan, Phys. Rept. 380, 235 (2003). 

E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006). 
E. V. Linder, Phys. Rev. D 80, 123528 (2009). 

P. Zhang, Phys. Rev. D 73, 123504 (2006). 

P. Zhang, M. Liguori, R. Bean, and S. Dodelson, Phys. Rev. Lett. 99, 141302 (2007). 

G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani, and S. Zerbini, Phys. Rev. 

D 77, 046009 (2008). 

K. Bamba, C. Q. Geng and C. C. Lee, JCAP 1008, 021 (2010) [arXiv: 1005.45741 [astro-ph.CO]] . 

Y. S. Song, W. Hu, and I. Sawicki, Phys. Rev. D 75, 044004 (2007). 

W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007). 

E. Komatsu et al. [WMAP Collaboration]. larXiv:1001.4538l [astro-ph.CO], 

S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 72, 123519 (2005). 

L. Perivolaropoulos, Phys. Rev. D 71, 063503 (2005). 

R. Amanullah et al, Astrophys. J. 716, 712 (2010). larXTv:1004,1711l [astro-ph.CO]. 

W. J. Percival et al, Mon. Not. Roy. Astron. Soc. 401, 2148 (2010). larXiv:0907. 16601 [astro- 

ph.CO]. 



D. J. Eisenstein and W. Hu, Astrophys. J. 496, 605 (1998), |arXiv:astro-ph/9709112 



E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009), 
larXiv:0803.0547l [astro-ph]. 

J. R. Bond, G. Efstathiou, and M. Tegmark, Mon. Not. Roy. Astron. Soc. 291, L33 (1997), 



arXiv:astro-ph/9702100 



W. Hu and N. Sugiyama, Astrophys. J. 471, 542 (1996), |arXiv:astro-ph/9510117 
A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002), [arXiv:astro-ph/0205436 
Y. Gong, X. ming Zhu, and Z.-H. Zhu. larXiv: 1008.50101 [astro-ph.CO]. 
K. Bamba, C. Q. Geng and C. C. Lee. larXiv:1007.0482l [astro-ph.CO], 



13 



[33] K. Glazebrook et al, ASP conference series 379, 72 (2007), |arXiv:astro-ph/0701876 



[34] D. Schlegel, M. White, and D. Eisenstein [with input from the SDSS-IH]. 1arXiv:0902.4680l 
[astro-ph.CO]. 

[35] G. J. Hill et al, ASP conference series 399, 115 f2008). larXTv:0806.0183l [astro-ph]. 
[36] European Space Agency Euclid Mission, |http: / / sci.esa.int /euclid/ 



[37] N. Gehrels, arXiv: 1008.4936 [astro-ph.CO]. See also http://jdem.gsfc.nasa.gov/ 
[38] D. J. Schlegel et al. larXiv:0904.0468l [astro-ph.CO]. 



[39] The Square Kilometre Array, http://www.skatelescope.org/ 



[40] J. A. Tyson, D. M. Wittman, J. F. Hennawi, and D. N. Spergel, Nucl. Phys. Proc. Suppl. 



124, 21 (2002), arXiv:astro-ph/0209632 J. A. Tyson [LSST Collaboration], Proc. SPIE Int. 



Soc. Opt. Eng. 4836, 10 (2002), arXiv:astro-ph/0302102 ; AIP Conf. Proc. 870, 44 (2006), 



arXiv:astro-ph/0609516 see also http://www.lsst.org/ 



[41] The Dark Energy Survey, http://www.darkenergysurvey.org/ 



11 



